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N2 pressure-broadened O 3 LINE WIDTHS AND 
STRENGTHS NEAR 1129.4 


By 

G.E. Copeland^, L.N. Majorana^ » C.N. Harvard^, 
and R.J. Sceinkamp** 


ABSTRACT 


A Beer 'a Law experiment waa performed with a tunable diode laaer to 
find the Nji preaaure-broadening characteriatica of a a ingle O3 abaorption 

line at 1129.426 cm“l (vl Band, J, Ka, Kc - 31, 1, 31 ♦ 32, 0, 32) for N2 
preaaurea from 10 to 100 torr <0^ preaaure - 3,16 torr). SO^ line poaitiona 
were uaed for wavelength calibration. Line ahapea were iteratively fitted 
to a Lorent* function. Reaulta were 6 (HWHM in MHz) - 47,44(±5.34) MHz 
1.730 (to, 088) MHz/torr *p(torr) with p ■ 0.9897. Thia intercept coraparea 
well with the Doppler Og - (^ broadened (at 3.16 torr) width of 44.52 Hz 
(ref. 1). Any difference i^j moat likely in laaer frequency uncertainty. 

Thia reaulta in a HWHM line width of 0.044 cm-^atm”^ at 760 torr and 285 K. 
The line atrengths integrated over Av "0.055 cm~l were found to be N2 prea- 
aure dependent; S(cm~2 torr"^ ) - i,334 x 10-*^p(torr)-^»3 with p ■ -0.9957. 


^Aaaociate Profeaaor, Department of Phyaica, Old Dominion Univeraity, 
Norfolk, Virginia 23508. 

^Graduate Reaearch Aaaiatant, Department of Phyaica, Old Dominion Univeraity 
Norfolk, Virginia 23508. 

^Reaearch Aaaiatant Profeaaor, Department of Phyaica, Old Dominion Univer- 
aity, Norfolk, Virginia 23508. 

^Undergraduate Reaearch Aaaiatant, Department of Phyaica, Old Dominion 
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INTRODUCTION 


u 


With the advent of semiconductor lasers and their recent application to 
spectroscopy (refs. 2-4) , there has been a revolution in the study of infra- 
red spectra. Use of the tunable diode laser (TDL) has perTsiilted hi^-reso- 
lution (<0.0001 cm-^ ) studies of several molecules of high current interest 

(CIO, HN(^ , SO2 i 0F2Cl2> etc.). This study is concerned with the experi- 
mental determination of the pressure-broadened characteristics of a single 

line near 1129.4 cm"', This line has been suggested as a possibi^ candi- 
date for remote sensing experiments using either infrared heterodyne spec- 
troscopy or differential lidar, since it is relatively well isolated from 
other atmospheric lines, has significant strength, and is relatively in- 
sensitive to temperature changes. A Beer's Law experiment is performed to 
determine the absorption coefficients and nitrogen (N2) broadening param- 
eters. 


APPARATUS AND PROCEDURE 

A simplified schematic of the apparatus is shown in figure 1. The 
radiation source is a PbSnTe TDL whose output is passed through collimating 
optics and a monochromator for mode selection. The output frcm the mono- 
chromator is divided by a beam splitter which directs 40 percent through a 

solid germanium etalon <>nto an LN2 -cooled HgCdTe detector. The remaining 60 
percent is split into 2 beams. One beam is directed through a reference gas 

cell to a second detector and the other through the 50-cm O3 absorption cell 
to a third detector. A more detailed description of the apparatus is to be 
found in the report by Major ana (ref. 5) . 

This arrrangement permits simultaneous recording of the O3 absorption 
spectra, the etalon fringes for relative frequency calibration, and the SO2 
reference gas spectrum for conversion to absolute frequency (ref. 4). The 
detector output, the TDL current, and S(^ and 0^ pressures were read by a 
12-bit A/D and transferred to a 7-track tape by a PDP8/e minicomputer for 
later analysis. O3 spectra were taken at 8 N2 pressures (10-100 torr) at 

285 K. SO2 et 1 torr was used as the reference spectrum since a detailed 
atlas exists (ref. 6) The 1129,418 cm"l line of 5O2 (vj ,V 2 ,V3 , j ,Ka,Kc) 
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[1,0,0,33,3,31 ^0,0,0,34,2,32] m g«lected «■ the standerd reference line 
becauae of its proximity to the 0 ^ line uwder atudy and the recent hetero- 
dyne confirmation of ita poaition at 1129f41835 

THE DATA AND ITS ANALYSIS 


The data were collected and proceaaed in the aame manner aa deacribed 
Majorana (ref, 5) . Strip chart recotdc were digitised and recorded on tape 
for later transfer to the ODU DEG System 10 computer, After input of data 
into the proper data format, each of the N 2 -broadened O 3 profiles was 
proceaaed using the program QUICK, One of the output products from program 
QUICK is a plot of the absorption coefficient as a function of wavenumber, 
as well as generation of data files containing the coefficient and the 
transmittance functions. 

For the foreign gas broadening study done in this research, the proper 
broadening function is a Lorenta function. Therefore, the Voigt function- 
fitting routine developed by Majorana and Copeland (ref, 5) cannot be used 
efficiently, A new program, BROAD (listed in the Appendix) was developed to 
fit in a nonlinear least squares sense a Lorentzian to the absorption coef- 
ficient: that is, a fit of the data is made to the function 


y 


A4 


A1 

(to-A2)2 


( 1 ) 


where A1 is the numerator of the Lorentzian, A2 is the center frequency 
(in cm“^), A3 is the half-width at half maximum (HWHM) of the profile, and 
A4 is the background. 


The program BROAD makes extensive use of the date-fitting routines 
developed by Bevington (ref. 8 ) and provides a file for plotting its fitted 
function together with the data for the user. The program begins by asking 
the user the name of the input data file. The program reports the mxnber of 
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data pointa, the date the data wap taken, the total (O 3 ♦ Ng) preaaure in 
torr/^ and the minlnium and maximum vnlut* of the absorption coefficient. The 
subroutine CHOP is called, and the user may reduce the size on high and low 
wavenumber ends of the spectrum to analyze. The total wavenumber interval 
across the data is reported, and the user may normalize the data if needed. 
The function AREA is called, and the integral of the absorption coefficient 
over wavenumber is calculated, i.e. , the strength in units of cm“^torr“^ . 

The user may select any of three different proc^:!!dures in which to wsight the 
data during the fitting procedure (MODE - +1,0,-1) which correspond to the 
Bevington weight convention (ref, 8 ). 

The program BROAD was designed to fit Lorentzians to either convex or 
concave data sets, i.e,, A4 near zero or AA maximum. The program auto- 
matically detects each case and, assuming only one spectral line, then makes 
initial guesses at the parameters Al, A2, A3, A4 and at their initial 

movements DELTAA(l-4). These are reported to the hser. The program uses 
the grid-search method, subroutine GRIDLS, and the DO loop 100 is entered. 
After each exit from GRIDLS, the reduced chi square statistic is compared 
with the previous value. If it has not changed by more than two percent, 
iteration is stopped and the program reports the final chi square, the num- 
ber of iterations, the final values of the parameters Al to A4, and 
estimates of their errors. Finally, the total pressure in torr and the 
half-width at half maximum in both cm**^ and MHz are reported. Ihe user may 
then select to see a recap of the data and the best fit and plot such as a 
graphics terminal. An example output of program BROAD is shown in the 
appendix. 

RESULTS 


These results are the first report to date of ozone line width broad- 
ening in the vl band using N 2 as the foreign gas. Table 1 lists the N 2 
pressure, the fitted half-widths, and the line strengths for the case where 
the initial ozone pressure was 3.16 torr. Figures 2 through 10 show the 
Lorentz fit to the absorption coefficient from the program BROAD. 
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The half'^widtht «• a function of preaaura ara shown in figure 11. Vfe 
find they can be represented by 


5(HWHM in MHa) ■ A + g p„ 


( 2 ) 


with 


A - 47.4 ± 6.3 MH* 
and 


B - 1.730 ± 0.01 MH*/torc 

with a linear correlation coefficient of 0.9897. Ihe intercept value 47.4 i 
5.3 MHz corresponds to the O 3 - O 3 broadened width at 3.16 torr 03 . This 
self-*broadening half-width compares well with the results of Majorana et al. 
(ref. 1), t’ho found a width of 44.52 MHz. The two values agree well within 
error bounds, Any significant difference could be traced to frequency in- 
stability in the diode laser. 

The strength of the absorption line is defined as the integral over 
wavenumber of the absorption coefficient; that is^ strength S is given by 

S<cm“2torr-M ■ k(v )dv (3) 

Generally, S is a constant for self -broadened data; however, we find S to 

be pressure dependent. Figure 12 shows the strength for this transition 

as a function of N 2 pressure. We find that the strength can be represented 

as a power function of the pressure, i.e., S ■ Ap§ , where A ■ 1.334 

"2 

X 10“^ and B ■ 1.314 ^ 0.049 with a correlation coefficient of -0.9957. 

Since -1.34 is essentially -4/3, then we find the strength varies as pres- 
sure of nitrogen raised to the -4/3 power. This curious result may be due to 
several effects: che^'ucal decomposition of O 3 due to collisions with N 2 or 

formation of w<C' 2 «. Van d* Waals molecules with N 2 , i.e. , O 3 + N 2 03 !^ and the 
resulting modification of residual ozone concentration. 
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APPENDIX 


PROGRAM BROAD 


Thif program ia designed perform the main data processing tasks as 
described in detail in the section "The Data and Its Analysis." This pro- 
gram is written in Fortran-10 and runs on a Dec System-10 time shared com- 
puter under the TopslO operating system. Extensive use is made of the 
following data handling routines selected from those of Bevington (ref. 3) 
and are not listed here; 


GRIDLS (X , Y .SIGMAY ,NPTS .NTERMS ,MODE ,A,DEI/TAA,SIGMAA , YPIT, 
CHISQR) 

INTEG ( X , Y ,NPTS , 1 , XI , X2 , SUM ) 

MAT1NV( ARRAY ,NTERMS,DET) 

FUNCTION FCHISQ( Y .SIGMAY ,NPTS »NFREE ,MODE . YFIT) 

FUNCTION AREA(X,Y,NPTS.NTERMS) 


LISTING 


PROGRAM BROAD 

WRITTEN BY 

DR. G. E. COPELAND 

DEPARTMENT OP PHYSICS AND GEOPHYSICAL SCIENCES 
OLD DOMINION UNIVERSITY 
NASA HETERODYNE GRANT 

DIMENSION YLABEL(i6),XLABEL(10) ,TITLE(35) ,SIGMAX(U00) , 
1SIGMAY(400) 

DIMENSION X ( 1*00 ) ,Y( 400 ) ,YPIT( 400 ) ,A(4 ) ,SIGMAA(4 ) 
DIMENSION DELTAA(4) 

DOUBLE PRECISION FILE 

SPEED OF LIGHT IN CM PER SECOND 

C-2.997925OEIO 

DATA FOR THE LABELS ON THE GRAPH 

DATA XLABEL/’W* ,’A' ,*V ,'E» ,’N’ ,'U» ,'M» ,*B' ,'E' ,'R»/ 


DATA YLABEL/’A»,'B'/S 

V,»BV 

DATA TITLE/»0S*3',» ' ,'L' ,♦!» .*11' ,»E» .» 'T» ,*1» , 'X% 

2,»T',*Z»,» '/P'.’I'/TV 
C 

DATA YIAB3/*R» ,»E» /S’ ,»X» ,»D’ ,*U’ ,’A* ,»L’ ,’S'/ 

TYPE 92 

92 FORMAT(’ LORENTZIAK FITTING FUNCTION’) 

WRlTE(5,XOOXO) 

ACCEPT X0020 ,FILE 

OPEN ( UNIT >2X .DEVICE- ^ DSK ’ .ACCESS- ’ SEQIN ’ .FILE-FILE) 
READ(2X.»)NPTS ,K2 ,K3 .K4 .P3 
IF(NPTS.GT.i400) TYPE 4X 

4X FORHAT(* TOO MANY DATA POINTS ENLARGE THE DIMENSIONS’) 
IF(NPT8.GT. 400) STOP 
TYPE 222.NPTS 

222 FORMAT {’ NUMBER OF DATA PAIRS -’ .l4) 

TYPE 45X.K2.K3.K4 

45X F0RMAT(XX.’ month -’.12.’ DAY -’.13.’ YEAR -’ .l4) 

C 

TYPE 45X9. P3 

45 X 9 FORMAT(’ PRESSURE IS ’.F7..*>' TORR’ ) 

C 

C READ THE DATA IN THE FILE INTO X AND Y ARRAYS 

C 

CALL READDT(X.Y.NPTS) 

C 

C 

DO 3X I-X.NPTS 
C TYPE 32 ,I.X(I).Y(I) 

32 FORMAT( XX,I3.XX,F8.3.XX.P8.3) 

3X CONTINUE 

C 

C 

CL0SE(UNIT-2X.DEVICE-’DSK’ .ACCESS^ ’SEQIN’ .FILE-PILE) 

C 

C 

CALL ROUTINE TO CHOP OUT XST AND LAST OP DATA TO GET RID OP OUTER LINES 

C 

C 

CALL CflOP(5PTS.X.Y) 

C 

C 

C 

C ASK IF DATA IS TO NORMALIZED 

C 

TYPE 77 

77 FORMAT/ DO YOU WANT THE DATA N0RMALIZED?<Y-X,N-0>’$) 

ACCEPT *,ANNOR 

IF ( ANNOR .EQ. O) GOTO 79 
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c 

C YES NORMALIZE IT 

C 

ANSW«AR£A(X , Y ,3 ) 

C 

C 

C SCALE THE Y ARRAY NOIf 

C 

DO 732 I»1»IIPTS 
732 Y(X)-Y(I)/A!ISW 

C 
C 

TYPE 765,ANSW 

765 FORMAT(» INTBQRAL OVER THE DATA IS -',1PE) 

79 OPEN ( UMIT-22 , DEVICE* ’ DSK • .ACCESS- » 8EQ0UT » .FILE-' PLTTEK, DAT » ) 

HSHET-1 
NGRP>2 

WRITE(22, 10190) N8HET 
C INPUT MODE fl.Q.-l. 

C 

C MODE-Kl IS THE INSTRUMENTAL 

C MODE - -1 IS STATISTICAL 

C MODE ——1 

C MODE -0 NO WEIGHTING 

C 

1*0 WRITE(5,10100) 

MODB--I 
ACCEPT •.MODE 

C CALCULATE MEAN IN XAND Y- 

YMEAN-O.O 
XMEAN-O.O 
DO 50 I-1,NPTS 
SIGMAX(l)-0.0 
XMEAN-X(I )-fXMEAN 
SIGMAY(l)-0.0 
YMEAN-YMEAN-I-Yd ) 

50 CONTINUE 

XMEAN-XMEAN/FLOAT(NPTS ) 

YMEAN-YMEAN /FLOAT (NPTS ) 

C CALCULATE 'SD* OP *Y' AND 'SD» OP *X» (SDY SDX) 

SDX-O.O “ 

SDY-O.O 

DO 70 I-l.NPTS 

f.DX- (X ( I ) -XMEAN ) ••2+SDX 

SDY»(Y(I )-YMEAN)»«2+SDY 

SIOMAX(I) - SQRT( (X(l)-XMEAN)»«2 ) 

SIGMAY(I)- SQRT((Y(I)-YMEAN)«»2) 

70 CONTINUE 

SDX-SQRT(SDX/PL0AT(NPTS-1 ) ) 

SDY-SQRT(SDY/PLOAT(NPTS-X ) ) 
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C WE HAVE SDY SDX TO PROVIDE START UP VALUES FOR 

C THE GRID SEAICK 

C 

YMAX-Y(l) 

DO 80 I»2,RPTS 

IF( YMAX .LT. Y{I) ) YMAX-Y(l) 

80 CORTIWJE 
NTERMS-4 
YMIN«Y(1) 

DO 81 I>2»NPTS 

IP( YMIM .or. Y(I) ) YMIH-Y(I) 

81 CONTINUE 
C 

C 

C GUESS ROUTINE 

C 

A(4)-( Y(1) +Y(NPTS) )*0.5 
C TYPE •,YNIN«YMEAN»YMAX,XMIN,XMEAN»XMAX 

A(2)-XMEAN 
C 

IF(Y(1) .LT. YMEAN ) GOTO 6ll 
C 

C CURVE IS CONCAVE 

ICON-1 
TYPE 877 

877 FORMAT {• CURVE IS CONCAVE') 

YMID«Y(l)-( Y(l)-YMIN)»0.5 

DO 992 I-l.NPTS 

IF( Y(I) .LT. YMID) XL»X(l) 

IF( Y(I) .LT. YMID) ILEFT*I 

IF( Y(I) .LT. YMID) GOTO 993 

992 CONTINUE 
STOP '992' 

993 DO 994 I-ILEFT,NPTS 

IF( Y(I).GT. YMID) XR*X(l) 

IF( Y(I) ,GT. YMID) GOTO 995 

994 CONTINUE 
STOP '994' 

811 CONTINUE 

C 

CURVE IS CONVEX 
TYPE 878 

878 FORMAT(' CURVE IS CONVEX') 

C 

YMID-Y ( 1 ) + ( YMAX-Y ( 1 ) ) «0 . 5 

DO 8121 I-1,NPTS 

IF( Y(I) .GT. YMID) XL«X(I) 

IF( Y(I) .GT. YMID) ILEFT*(I 
IF( Y(l) .GT. YMID) GOTO 8131 
8121 CONTINUE 

STOP '8121* 


DO 8l4 I«»ILEFT,NPrS 

IP{ Y(I) .LT. YMID) XR-X(I) 

IP( Y(I) .LT. YMID) GOTO 995 

CONTINUE 

STOP »8l4' 

A(3>XR-XL 

A3 IS FWHM 

A(l)-(0.25» A(3)»»2+(XMEAN-A(2))»»2)«( YMAX-A(U) +0.5) 
IF( Y(1).GT. YMIN) A(1)«-A(1) 

GOT A1,A2,A3,A4 


TYPE *,A 

NOW HAVE FINISHED THE 

GENERATION OF GUESSES ON A1,A2,A3,A4 FOR A LORENTZIAN FUNCTION 
Y-A4 + Al/( ( X-A2)«»»2+( A3/2)««2) 


WRITE(5*10120) 

WRITE ( 5 , 10110 ) ( I , A( I ) , 1*1 , NTERMS ) 

C FINISHED STARTING VALUES 

C NOW FIND THE INCREMENTS ON A(i) ,A(2) ,A(3) 

C SET STEP SIZE IN MEAN * 6SDX/20 

DELTAA(3)*A(3)/10.. 

C SET STEP SIZE ON COEFFICIENT A1 

DELTAA(l)-A(l,)/10. 

C SET STS? SIZE ON COEFFICIENT A2 

DELTAAC2 )■ (XMEAN+SDX ) /lO .0 

C GIVE GUESS ON DELTAA(4) BACKGOUND 

DELTAA(4)*A(4)/10.0 

C 

C ENTER LOOP 

C ITERATE 20 TIMES AT MOST 

kNT»0 

DO 100 K*l,15 

CALL GRIDLS ( X , Y , SIGMAY , NPTS , NTERMS , MODE , A ,DELTAA , 

1 SIGMAA.YFIT.CHISQR) 

TYPE 88T,K.CHISQR 

IF(K .GT. 1) GOTO 90 

SAVE*CHISQR 

SAVIT*SAVE 

KNT*1 

GOTO 100 

90 XCHI»100.0*(CHISQR-SAVE)/CHISQR 
IF( XCHI .GT. 0.0 ) GOTO 110 

c 

CRITERION OF CONVERGENCE IS CHISQR DOES NOT CHANGE BY 1 PER CENT 




c 

C OR IP CHISQR IS LESS THAN 0.02 OP THE 1ST CHISQR POUND 

C 

C BEVIPOTON PAGE 212 

C 

IP( ABS(XCHI) .LT. 0.5 ) GOTO 110 
IP( CHISQk'i .LT. 0.020«SAV1T) GOTO 110 
SAVE - CHISQR 
KNT-KNT+1 
100 CONTINUE 

110 KNT-KNT+1 

WRITE (5, 10130) KNT 

C RETURNED VALUES ARE A, SIGMAA.YPIT, CHISQR 

WRITE(5,10140) 

DO 120 K.NTERMS 

WRITE( 5 ,10150 ) ,I ,A( I ) ,SIGMAA( I ) 

120 CONTINUE 

WRITE (5 , 10160 ) , CHISQR 
WRITE(5. 10161 ),P3 

10161 PORMAT(» TOTAL PRESSURE IS ■»,P6.2,* TORR') 

A3K-A(3)».5 
A3MHZ*A3K*C/1.E+6 
TYPE 345 ,A3K,A3MHZ 

345 FORMAT ( IX, 'HWHM ■•,P9.6,» CM-1 OR ',F12.4' MHZ’) 

130 WRITE (5, 10170 ) 

C ANS1»0 

ACCEPT »',ANS1 
IP(ANS1.EQ.1)G0 TO I 50 
l40 WRITE(5,10180) 

ACCEPT »,ANS2 

IF(ANS2.EQ.O) STOP 

WRITE( 22, 10190 ) NGRP 

WRITE(22, 10190 ),NPTS 

WRITE( 22 . 10080 ) ( YLABEL (I ) ,I» 1 , 16 ) 

WRITE<22 , 10080 ) (XLABEL (l ) ,1*1 ,10 ) 

WRITE ( 22 , 10200 ) ( TITLE ( I ) ,1*1 , 35 ) 

WRITE( 22 , 3451 ) ,A3K , A3MHZ , P3 

3451 FORMAT ( IX, 'HWHM ',F9.6,» CM-1 ',P7.2' MHZ P» »F6.2* T0RR») 

C OUTPUT TO TEKTRONIX PLOTTING FILE 

C 

DO 812 I*1,NPTS 

WRITER 22, 10220 ) X(I),Y(I) 

812 CONTINUE 
C 

C WRITE FITTED DATA AND GENERATE NEW X VALUES 

C 

C 

C WRITE OUT THESE X AND YFIT 

DO 813 I»1,NPTS 

813 WRITE(22, 10220) X(l),YFIT(l) 

99 CONTINUE 

I 
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CLOSE ( UNIT-22 , DEVICE- ’ DSK ’ .ACCESS- » SEQOUT ' .PILE- » PLTTEK. DAT ’ ) 

STOP 

C OUTPUT A DATA MATRIX IN TERMINAL 

150 CONTINUE 

Q-0.0 

C 

C FIND AVG OF ABS{% ERROR) 

C 

WRITE(5,10240) 

DO 88 I-l.NPTS.lO 

z-( (x(i)-YPiT(i))/mT(i) )noo.o 
Q-Q-»-ABS(Z) 

88 WRITB(5..10250) X(l).Y(l) .YPIT(I).Z 

Q-10.»Q/PL0AT(NPTS) 

WRITE (5. 10230 )Q 
WRITE (5. 10240) 

GO TO l40 

887 F0RMAT(» FINISHED ITERATION ’.12.’ CHI SQ -».1PE) 

10010 FORMAT (IX. 'WHAT IS INPUT PILE NAME<CONTAINS OUTPUT PROM 
QUICK>?') 

10020 FORMAT (AlO) 

10060 PORMAT(I) 

10080 FORMAT (40A1) 

10090 FORMAT (1X,40A1) 

lOlOO PORMATUX* 'WHICH MODE DO YOU WISH TO PIT DATA*,/. 

1 IX , ' ■►1-INSTROMENTAL PITS- 1 /VARIANCE ' , / . 

2 lX.'O- NO WGHT’ ,/, 

3 IX.* -1- STATISTICAL-l/Y » *,/,lX,*?*$) 

10110 P0RMAT(5X,* A(*,II.* ) »*,1PE) 

10120 P0RMAT(10X,* START UP VALUES*) 

10130 FORMAT ( 1 X , * THERE WERE * , 13 . ' ITERATIONS * ) 

10140 P0RMAT<1X. 'USING Y-a4+A1/( (X-A2)**»2 + 0.25'*A3 ""S) 

I'./.IX. 'YOU OBTAIN: *) 

10150 F0RMAT(5X,*A(’.I1,*) -* ,1PE, *•»•/- ' .IPE) 

10160 FORMAT (5X, 'CHI SQUARED - *,1PE) 

10170 P0RMAT(1X,/,1X,'D0 YOU WANT A DATA RECAP? YES«l,N0-07'$) 

10180 P0RMAT(2X,'D0 YOU WANT TO PLOT DATA ON GRAPHICS ?( YES-1, NO-0 ) '$) 

10190 FORMAT (2X,1) 

10200 FORMAT (4oA1) 

10210 PORMAT(5X,*Al-4-*,4(lX.lPE9.2) ) 

10220 FORMAT (2X, IPE, 2X, IPE) 

10230 FORMAT (IX, 'MEAN OP THE ABSOLUTE VALUE OP THE % ERROR -' ,P) 

10240 F0RMAT(2X,/,/,8X,»X-DATA' ,8X,*Y«DATA' ,8X, 

I'YPIT'.aX.'PER CENT DIFP') 

10250 FORMAT ( IX , 3 (IPE , 3X ) , OPPIO . 2 ) 

END 

C 

C— 

FUNCTION FUHCTN(X,I,A) 

C 

C 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

10 

30 
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PURE LORENTZIAN FUNCTION 
Y»A4+(A1)/( (X-A2)«»2 +(A3/2)**2 ) 


PLUS A BACKGROUND 

DIMENSION X(l),A(l) 

XI*X(I) 

Z»XI-A(2) 

Z2*Z*Z 

PUNCTN«A(1)/(Z2+(A(3)*A(3)»0.25)) 

FUNCTN*FUNCTN+A{4) 

RETURN 

END 


SUBROUTINE READDT(X,Y,N) 

DIMENSION X(l),Y(l) 

READ (21 ,211 )KK ,N ,K2 ,K3 ,K4 ,K10 ,TAPENO , P3 
F0RMAT{1X,I2,1X,I4,1X,3(I2,1X),I1,1X,A5,1X,P) 

READ( 21 ,EJn)»22 . * ) ( ( X ( I ) , Y ( I ) ,I*J , J+4 ) , J»1 ,N , 5 ) 

F0RMAT(10(1X,1PE8.2) ) 

WRITE(5,334)( (X(I),Y(I),I-J.J+4),J*1.N,5) 

RETURN 
CONTINUE 
TYPE 33 

FORMAT (' RAN OUT OF DATA IN READDT RECAP FOLLOWS') 

TYPE 34 *N 

FORMATC number of data POINTS ■',!) 

WRITE(5,*)( (X(i),Y(I),W,J+4),J-1*N,5) 

RETTURN 13 


END 

C 

C — 

C 

SUBROUTINE CHOP(N,X,T) 

DIMENSION X(l),Y(l) 

DIMENSION XSAV(1*00),XSAV(1»00) 

C 

C CHOP OUT PARTS OP THE ARRAY AND REPLACE IT IN THE ORDINAL ARRAYS 

C 

NSAV-N 

DO 10 I-l.NSAV 
XSAV(I)-Xd) 

10 YSAV(I)-Yd) 

YMIN«Y(l) 

DO 20 I«2,NSAV 

IP(YMIN .GT.Y(I)) IMIN-I 

IP(YMIN .GT.Y(I)) YMIN«Y(I) 

20 CONTINUE 

TYPE 100 ,YMIN,IMIN 

100 F0RMAT(» MINIMUM VALUE OF Y *'P,' AT POINT ',13) 

YMAX«Y(l) 

DO 201 I-2.NSAV 

IP( Y(I) .GT. YMAX) IMAX»I 

IF( Y(I) .GT. YMAX) YMAX-Y(I) 

201 CONTINUE 

TYPE 2011, YMAX, IMAX 

2011 FORMAT (' MAXIMUM VALUE OP Y »'F,' AT POINT »,I3) 

109 TYPE 110 

110 FORMAT(' DO YOU WANT TO CHOP PARTS OF THE ARRAY OFF7<Y»1,NO»0>' ) 
ACCEPT »,ANS 

IF(ANS.EQ.O) RETURN 
IF(ANS.EQ.l) GOTO 150 
nOTO 109 
150 CONTINUE 

TYPE l60 

160 FORMAT (' INPUT OF LEFT DATA POINT TO START NEW ARRAY? '$) 

ACCEPT »,IST 

TYPE l6l 

161 FORMATC input OF DATA POINT* ON RIGHT TO STOP THE ARRAY? »$) 
ACCEPT »,ISP 

IF( ISP. LT. 1ST) GOTO 109 
C 

C REPORT THE DIFFERENCE IN FREQUENCY BISTWEEN THE 2 LIMITS 

C 

DIFF«X(ISP)-.X(IST) 

TYPE 1601,DIFF 

1601 FORMAT (• WAVENUMBER CHANGE ACROSS THIS SAMPLE «',F8.7,‘ CM-1') 

DO 30 I*IST,ISP 
Xd-IST+l)»XSAV(l) 
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Y(I-IST+1)»YSAV(I) 

NilSP-IST 

RETURN 

END 


' RUN BROAD 


^ « 

LORENTZ IAN FITTING FUNCTION 

WHAT- 18 INPUT FILE NAHE<CONTAINS OUTPUT FROH EARL>T 

FILS ‘ * 

NUHBER OF DATA PAIRS -200 
MONTH -12 DAY ■ 19 YEAR -79 
PRESSURE IS 53*59 TORR 

MINIMUM VALUE OF Y » -0*0000024 AT POINT 180 

MAXIMUM VALUE OF Y - 0*0000723 AT POINT 93 

DO YOU WANT TO CHOP PARTS OF THE ARRAY OFFTCY-A rN0*0> 
1 

INPUT ♦ OF LEFT DATA POINT TO START NEW ARRAY? 

20 ' 
INPUT f OF DATA POINT TIN RIGHT TO STOP THE ARRAY? 

180 

WAVENUMBER CHANGE ACROSS THIS SAMPLE -*0484009 CM-1 
DO YOU WANT THE DATA N0RMALIZED?<Y»1 rN-0> 

1 * . 

INTEGRAL OVER THE DATA IS - 8* 1862793E-07 

.WHICH MODE DO YOU WISH TO FIT DATA 
+1-INSTRUMENTAL FITS- 1/VARIANCE 
0- NO WGHT 

-1- STATISTICAL-l/Y - 


0 

CURVE IS CONVEX 

START UP V A 
A<1 > « -l*5973438E-03 
A<2 ) - l*1294275E+03 

A<3 > - 8*4991455E-03 

A(4 > - 3*6646686E-01 

FINISHED' ITERATION * 1 CHI SG 
FINISHED ITERATION # 2 CHI SG 
FINISHED ITERATION * 3 CHI SG 
FINISHED ITERATION ♦ 4 CHI SO 
FINISHED ITERATION * 5 CHI SG 
THERE WERE 5 ITERATIONS 
USING Y»A4+Al/< <X“A2)*#2 + Q. 
YOU obtain: 


LUES 


- l*2879837E+<)2 

- l*275593AE+02 

- l*2667750E+02 

* l*2597696E+02 

- l*2537410E+02 

25JKA3 ##2) * 


Ad) - l*6568207E-03+/- -4 ♦ 1651741E-06 

A(2> - -l*1294275E+03+/- 3*5054553E-03 

A (3) « 9*0509B54E-03+/- 1 .4426374E-05 

A<4) » -3*5671993E-01+/- 7 *905A752E-02 

CHI SQUARED - 1 *2S37410E+02 ‘ 

TOTAL PRESSURE IS - 53*59 TORR . 

HWHM = 0*004525 CM-1 OR 135*6709 MHZ 



* 


■ , *■ 

DO YOU U^NT A DATA RECAPT YCS-l»NO-OT 
0 

DO YOU WANT TO PLOT DATA ON GRAPHICS ?( YES-1 rN0-0> 

0 

STOP » . 

END OF EXECUTION 

CPU TIKEr 1*46 ELAPSED TIHEt n3.32 
• EXIT 
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Table I* Half-widcha (HWHM) and strengtha va. Ng preaaure, 

for O 3 broadened with O 3 preaaure ■ 3,16 torr at 
285 K. 


SI^ Preaaure 
(torir) 

Half 

(mKayaera) 

-widtha 

(MHc) 

Strengtha* 
0 j 

cm“^torr ‘ 

0.0 

1.615 

48.43 

4.602C-5) 

10.15 

1.929 

57.83 

6.738C-5) 

15.10 

2.290 

68.66 

4.043C-6) 

19.95 

2.784 

83,47 

2.385(-6) 

30.30 

3.819 

114.5 

1.435(-6) 

50.43 

4.491 

134,6 

7.892(-7) 

60, 4S 

4.945 

148.3 

5.945(-7) 

90,80 

7.600 

227.8 

3,853(-7) 

100.6 

6.834 

205.0 

3.438(-7) 


^integration over 55 mKayaera 
°Numbera in v ) are powers of 10 
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TDL SPECTROSCOPY FACILITY 

I. R. ABSORPTION CELL HgCdTe DETECTOR 



Figure 1. Optical setup for simultaneous measurement of reference gas 
spectnnn, test gas spectrum and etalon tuning curve. 
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Figure 6 
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Absorption coefficient vs. wavenumber at pressure - 33.46 torr. 
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Figure 7. Absorption coefficient vs. waveniraber at pressure - 53.59 torr 
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Figure 8, Absorption coefficient va. wavenumber at pressure ■ 63.64 torr 
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Figure 9s Absorption coefficient vs. wavenumber at pressure - 93.80 torr 






